Multiple-spin coherence transfer in linear Ising spin chains and beyond: 
numerically-optimized pulses and experiments 
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We study multiple-spin coherence transfers in linear Ising spin chains with nearest neighbor cou- 
plings. These constitute a model for efficient information transfers in future quantum computing 
devices and for many multi-dimensional experiments for the assignment of complex spectra in nu- 
clear magnetic resonance spectroscopy. We complement prior analytic techniques for multiple-spin 
coherence transfers with a systematic numerical study where we obtain strong evidence that a cer- 
tain analytically-motivated family of restricted controls is sufficient for time-optimality. In the case 
of a linear three-spin system, additional evidence suggests that prior analytic pulse sequences using 
this family of restricted controls are time-optimal even for arbitrary local controls. In addition, 
we compare the pulse sequences for linear Ising spin chains to pulse sequences for more realistic 
spin systems with additional long-range couplings between non-adjacent spins. We experimentally 
implement the derived pulse sequences in three and four spin systems and demonstrate that they are 
applicable in realistic settings under relaxation and experimental imperfections — in particular — by 
deriving broadband pulse sequences which are robust with respect to frequency offsets. 



PACS numbers: 82.56.-b, 03.67.Ac, 02.30.Yy 



I. INTRODUCTION 

The control of spin dynamics in chains of coupled 
spins- 1/2 is a topic of both theoretical and practical in- 
terest |l|-[l2l|. On the one hand, the use of spin chains 
is considered for the efficient transfer of information in 
future quantum computing devices [l3l.[l^. On the other 
hand, coherence transfer between remote spins is the ba- 
sis of many multi-dimensional experiments for the assign- 
ment of complex spectra [Tsl - flTj in nuclear magnetic res- 
onance (NMR) spectroscopy. In addition to linear spin 
chains with only nearest neighbor couplings, in realistic 
settings also long-range couplings between non-adjacent 
spins must be considered. For example in ^"^C and ^^N 
labeled proteins, the nuclei in the protein backbone form 
a chain of coupled spins- 1/2 with dominant next neigh- 
bor (single bond) couplings and smaller and 
couplings (via two or three chemical bonds) between non- 
adjacent spins in the chain [17]. 

Here we focus on the efficient creation of multi-spin 
operators from a single-spin operator in a spin chain, such 
as the creation of multiple-spin order from polarization 
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of the first spin 

hz ^ 2" ^Ilzl2z • • • I{n-l)zlnz- 

The transfer shown in Eq. ^ is just a prototype example 
of a general transfer of the form 



(1) 



hs 



(2) 



where 5, efc G {x,y,z} for k = l,...,n. Note that the 
transformations in Eqs. ((!]) and ^ are identical up to lo- 
cal spin rotations. Hence, in the limit where the time for 
selective rotations of individual spins is negligible (com- 
pared to 1/ {2Jmax), where Jmax is the largest spin-spin 
coupling constant in the chain), the transformations in 
Eqs. Ill) and ([2]) can be achieved in the same amount of 
time. (This situation is typical for heteronuclear NMR 
experiments in liquid state, where the control amplitudes 
for single spin operators are orders of magnitude larger 
than the largest coupling constants.) In Eq. ([2]), the ini- 
tial single-spin state is not limited to longitudinal mag- 
netization (polarization Ji^) but may also be transverse 
magnetization (in-phase coherence hx or Iiy) [13, fisj . 
Examples of multi-spin target operators in Eq. ([2]) con- 
taining one or several transverse operators include states 
of the form 2"^^^ Ii^hz ■ ■ ■ I(n-i)zlnx (corresponding to 
anti-phase coherence of spin n with respect to spins 1 
to n - 1), and 2'^~^Iixl2x ■ ■ ■ I{n-i)xlnx (corresponding 
to multi-quantum coherence), which are relevant in so- 
called "out and back" transfer schemes [13, [H, HO] and 
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in the creation of multiple-quantum coherence [18|, |2l[ , 
respectively. 

We consider the case of Ising-type spin chains [13, [l^ , 
corresponding to the weak coupling limit [l8| . In this 
limit, which is an excellent approximation for hetero- 
nuclear NMR experiments, the coupling Hamiltonian for 
a pair of spins k and I has the form 

where Jki is the coupling constant in units of Hertz 
(Hz). In conventional experiments, the standard meth- 
ods to achieve transfer in Eq. are based on COSY- 
or RELAY- type transfer steps [ITI ITsj . which are real- 
ized in hetero-nuclear experiments by a series of INEPT 
building blocks [13] (see Sec. |TT|. The transfer time is 
determined by the size of the coupling constants J^i in 
a given spin system. For example, in a linear spin chain 
with only next-neighbor couplings, the total duration is 
given by 



Assuming a restricted pulse structure (see Se es. IIII Al 



We are interested in finding the shortest possible time to 
achieve the transfer in Eq. ([2]) or conversely, the max- 
imum transfer amplitude for any given time, which re- 
mains an open question up to now. 

For relatively simple spin systems, consisting of up 
to three spiiis, ^time-optimal j25i431| and relaxation- 
optimized [32h37I | pulse sequences have been recently 
found analytically based on methods of optimal control 
theory [UlUI, establishing rigorous physical limits of 
minimal transfer times or minimal relaxation losses, re- 
spectively. In addition to powerful analytical tools, op- 
timal control theory also provides efficient numerical al- 
gorithms for the optimization of pulse sequences, such 
as the gradient ascent pulse engineering (GRAPE) algo- 
rithm, exploiting the known equation of motion for the 
spin system j4^3l - |46j . With this algorithm it is possible to 
optimize tens of thousands of pulse sequence parameters 
and the resulting pulse sequences are not limited to pre- 
viously known transfer schemes. However, in contrast to 
analytical methods proving global optimality of a given 
pulse sequence, there is no guarantee that numerical op- 
timal control algorithms like GRAPE will converge to 
the global optimum [IJ- Never-the-less, in cases where 
the theoretical limits are known, the GRAPE algorithm 
closely approached these limits [4^ . This motivated 
its use also in cases for which analytical results on the 
global optimum are presently unknown in order to ex- 
plore the physical limits of the maximum possible trans- 
fer efficiency as a function of transfer time, resul ting in so- 
called TOP (time-optimal pulse) curves [41|,|49|-l53|[. Fur- 
thermore, additional effects such as relaxation |54l. [ssj . 
radiation damping jsoj , and experimental constraints and 
imperfections — such as limited control amplitudes and 
control field inhomogcneitics [57l - [59j — can be taken into 
account to find highly robust pulses suitable for practical 
applications under realistic conditions. 



and IIV A)) analytic pulses were derived in Refs. and 
[l^ , respectively, for the case of equal and unequal cou- 
plings. This results in significantly shorter transfer times 
compared to conventional approaches, however it was not 
clear how closely the performance of the derived pulse se- 
quences converges to the time-optimal performance. 

In this work, we summarize the analytic approach of 
Refs. [H [il (see Sections HITXl and and explore 

its time-optimality by conducting a systematic numeri- 
cal study of the considered coherence transfer (see Sec- 
tions HTlB] and |IVB|). Focusing on the case of linear Ising 
spin chains with three and four qubits we compare the 
duration of pulse sequences for arbitrary pulse structures 
with the restricted pulse structure motivated by the an- 
alytic pulses. We also numerically analyze linear Ising 
spin chains for up to six spins (sec Section IVIII) . In ad- 
dition, our numerical approach makes it also possible to 
investigate more realistic spins systems with more gen- 
eral coupling topologies (see Section |V| . 

In Sec. I VII we show how to make the pulse sequences 
robust with respect to off-resonance effects using the 
DANTE (Delays Alternating with Nutations for Tailored 
Excitation) approach [2^ [6^ |6l|. Finally, we present 
experimental results for model spin chains consisting of 
three and four hetero-nuclear spins-1/2, demonstrating 
good performance of the new sequences under experimen- 
tal conditions and comparing the results to conventional 
pulse sequences. 



II. COHERENCE TRANSFER IN LINEAR 
ISING SPIN CHAINS 

Throughout this work we mostly consider linear Ising 
spin chains which have only direct couplings between 
neighboring spins [13, [1^ . (Later we will also allow ad- 
ditional couplings between non-neighboring spins.) As- 
sume that a chain of n spins is placed in a static external 
magnetic field along the z-dircction and that neighboring 
spins are coupled by an Ising interaction where the cou- 
pling strengths Ji^i+i are fixed but may depend on the 
position 1 < £ < n—1 in the chain. Without any control, 
the system evolves freely under its drift Hamiltonian 



Hd = 2 TT ^ Jui + ilizl(i+i)z- 
1=1 

The drift Hamiltonian is given in a suitably chosen mul- 
tiple rotating frame, which rotates simultaneously at the 
resonance frequency of each spin. We use the product- 
operator basis Iiv= 0, where aj=v for j=^ and aj=0 
otherwise (see Ref. [i3)- The matrices Ix'=[W) /'i,, 
Iy:= ( ° ) /2, and /^:= ( J i'j ) /2 are the Pauli spin ma- 
trices and ^o:=(oi) is the (2 x 2)-dimensional identity 
matrix. In addition to the free evolution, we assume that 
individual spins can be selectively excited using radio- 
frequency (rf) pulses, which is the case if the Larmor 
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FIG. 1. A linear three-spin chain has only direct couplings 
Ji2 and J23 between neighboring spins. 



frequencies of the spins are well separated as compared 
to the coupling strengths Ji.i+i- Thus controls on indi- 
vidual spins can be applied on a much faster time scale 
as compared to the free evolution w.r.t. the drift Hamil- 
tonian. 

We derive explicit controls for the amplitude and phase 
of the external rf fields implementing a unitary evo- 
lution which transforms an initial polarization Ii^ on 
the first spin to a multiple-spin state 2"~-'^(]^"^j^ Iiy)Inz 
while minimizing the pulse duration tp. In the following, 
we often compare control pulses with the conventional 
strategy, which consists of n—1 steps of free evolution 
(1 < TO < n—1) 

r)m-lr r r nm j j j 

^ ^ly ' ' ' ^m—l,y^mx ^ ^ ^ly ' ' ' ^my-'-ra+l^z 

where each individual step — besides the final one — is fol- 
lowed by one hard -j-pulsc on the (m+l)th spin along 
the y-dircction. As each period of free evolution is of 
length l/{2Je^e+i) where Jej+i is given in Hz, the total 
evolution time is given by tp = l/(2t/£,£+i)- 




2 4 6 8 9.8 5 10 15.5 



t (ms) / (ms) 

FIG. 2. Analytic pulses for linear three-spin chains are given 
in the cases of (a) k = J23/J12 = 88.05/88.05 = 1 and (b) 
k = 1.59 ^ J23/J12 = 73.1/46. 




FIG. 3. The schematic coupling topologies of (a) ethanamide 
and (b) diethyl-(dimethylcarbonyl)fluoromethylphosphonate 
(see [12, H^l) result in experimental three-spin systems with 
coupling ratios (a) = 1 = 88.05/88.05 and (b) k = 1.59 « 
73.1/46. Larger couplings are shown as solid-black lines, and 
smaller couplings are shown as dashed-black lines. Decoupled 
spins are given in gray color. 



III. LINEAR THREE-SPIN CHAINS: 
ANALYTIC AND NUMERICAL APPROACHES 

A. Analytic approach 

In this section, we consider the model of Sect. [ll] in the 
case of linear three-spin chains (see Fig. [T]). In the most 
general case, one could allow independent controls on 
each of the three spins along both the x- and y-direction. 
But in order to simplify the control problem we allow 
only one control on the second spin along the y-direction. 
This might not lead to time-optimal controls. But even 
using this restricted model, controls which are shorter as 
compared to the conventional strategy were obtained in 
Ref. [l3| (see also [HI). In the following, we summarize 
the analytical approach of Ref. [l^l . 

Starting from an initial state Iix and using only one 
control on the second spin along the y-direction, we can 
analyze the control problem on the subspace spanned 
by the operators Iix, iliyhz, 2Iiyl2x, and iliyhylsz 
as compared to the full 63-dimensional space of opera- 
tors. Using the notation (O) := Tr(Op) for the expec- 
tation value and Tr for the trace, we denote the cor- 
responding expectation values by xi = xi{t) ~ {hx), 
X2 = X2{t) = {2Iiyl2z), X3 = X3{t) = {2Iiyl2x), and 
X4 = X4{t) = {^Iiyhyhz) ■ We obtain the differential 
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where u — u{t) denotes the amplitude of the control on 
the second spin along the y-direction and k ~ J23/Ji2- 
Using the coordinates (xi, a;2, X3, 2:4)"^ we aim to time- 
efficiently transfer (1, 0, 0, 0)^ to (0, 0, 0, 1)^. 

Now, we change from the coordinates {xi,X2,X3,X4)'^ 
to the coordinates 



(fl,f2,?-3)^ = {Xl,\jxl +xl,Xif 

on the sphere where 9 = 9{t) is given by tan0 = xz/x2- 
This transforms Eq. ^ to 

fi\ / -cos6l \ /ri\ 
7-2 — ^ cosO — fcsin6' M'2 ■ 
fzj \ fcsinfl j \r3j 

In the new coordinates, we want to time-efficiently trans- 
fer (1,0,0)^ to (0,0,1)^. 

In order to find the time-optimal controls, Euler- 
Lagrange equations were set up and solved in Ref. [l^ 



4 



TABLE I. We compare the duration tp, the logarithmic fi- 
delity F, and the shape of numerically-optimized pulses for 
a linear three-spin chain with coupling ratios (a) fc = 1 and 
(b) k = 1.59. The number of controls #u is given in the first 
column. In the third column we present the corresponding 
logarithmic TOP curves. The second column shows an exam- 
ple of a shaped pulse, whose position is denoted with an x in 
the logarithmic TOP curve. The rf control on the middle spin 
along the y-axis is plotted using a solid-black line. Other rf 
controls are plotted using dashed or solid-gray lines. 



pulse shape 



logarithmic TOP curve 




(b) k = 1.59 
^60 




leading to the differential equation 



1 



l26I 



(4) 



for the variable 6. The differential Eq. Q ean be nu- 
merieally integrated if the initial values 9{0) and ^(0) 
are known. Using the results of Ref [l^ one can de- 
termine conditions on the initial values: In the case of 
(ri(0),r2(0),r3(0))^ = (1,0,0)^ one can deduce that 
6(0) = 0, but 9(0) is undetermined. In Ref. [13] combi- 
nations of one-dimensional searches were used to deter- 
mine the optimal 0opt(O the time-optimized control 
as Wopt(i) = Ji29opt(t). Examples for the corresponding 
(semi-) analytic pulses are shown in Fig.O The values are 
motivated by the experimental systems given in Fig. [3] 



TABLE II. For coherence transfers in linear three-spin chains 
(fc = 1 and k = 1.59), we give the numerically-optimized 
times tp and the fidelities F in the case of one, two, and six 
rf controls (see text). The duration tp is independent of the 
number of controls which suggests that only one rf-control on 
the middle spin is sufficient for the time-optimal coherence 
transfer. 



k 


#u 


tp (s) 


F 


k 


#w 


tp (s) 


F 


1 


1 


0.0098 


0.9999 


1.59 


1 


0.0155 


0.9999 


1 


2 


0.0098 


0.9999 


1.59 


2 


0.0155 


0.9999 


1 


6 


0.0098 


0.9999 


1.59 


6 


0.0155 


0.9999 



steer an initial state (or unitary transformation) to a final 
state (or unitary transformation) while minimizing (e.g.) 
the duration of the pulse. Both the amplitude and the 
phase of the resulting pulse can be smoothly- or noisy- 
shaped depending on e.g. the initial pulse or bounds on 
the control strength (see, e.g. [ssj). 

We treat three different levels of rf controls: First, we 
use only one rf control operating on the second spin along 
the j/-direction. Second, we use two different rf controls 
operating on the second spin along both the x- and y- 
direction. Third, we use a total of six rf controls oper- 
ating on each of the three spins along both the x- and 
y-direction. We remark that employing rf controls on 
one spin along both the x- and y-direction gives com- 
plete (local) control on that spin. Let k denote the ratio 
between the couplings J23 and J12. We determine the 
numerically-optimized pulses and plot the logarithmic fi- 
delity F vs. the duration tp of different shaped pulses 
for the coupling ratios fc = 1 and fc — 1.59 which are 
motivated by the experimental scenarios of Fig. |3l The 
numerical results are given in Table ID We show exam- 
ples of shaped pulses with short duration tp and high 
fidelity F > 0.9999. In addition, we present logarithmic 
time-optimal (TOP) curves where we plot the logarith- 
mic transfer efficiency (i.e. log(l— F) where F is the fi- 
delity) versus the optimal transfer time. Comparing the 
different cases suggests that only one rf control on the sec- 
ond spin is sufficient for a time-optimal pulse. For high 
fidelities {F > 0.9999), the durations of the analytic and 
numerically-optimized pulses are identical (in the given 
accuracy) while the pulse forms differ. In Table [TTl we 
compare the duration of pulses on linear three-spin sys- 
tems for different values of fc. 



B. Numerical approach 

We numerically optimize pulse shapes by employing 
the GRAPE algorithm [l^ which was developed by em- 
ploying principles of optimal control theory. Using a 
gradient-based optimization we obtain rf controls which 



Conjecture 1 Consider a linear three-spin chain with 
local controls on each spin. One can time- optimally 
transfer coherence from Iix to '^Iiyl2yl5z using only one 
control on the second spin along the y-direction. In ad- 
dition, the analytical pulses of Refs. fJH . [7^ / are time- 
optimal in the case of linear three-spin chains even if 
one allows arbitrary local controls. 
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FIG. 4. A linear four-spin chain has only direct couplings 
Ji2, J23, and J34 between neighboring spins. We split the 
corresponding four-spin chain control problem into two sub- 
problems for three-spin chains. 
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FIG. 5. Analytic pulses for linear four-spin chains are given in 
the cases of (a) fci = fc2 = 1 as well as (b) k\ — 2.38 and fc2 — 
0.94. The pulses on the second and third spin along the y- 
direction are given, respectively, as solid and dashed line. The 
corresponding two hard pulses on the second and third spin 
are depicted by a vertical line where the flip angles are given 
above. The hard pulses in the left figure can be implemented 
by applying a pulse of 5000 Hz for 17.40 microseconds. The 
hard pulses in the right figure can be implemented by applying 
a pulse of 5000 Hz for 11.21 microseconds. 



IV. LINEAR FOUR-SPIN CHAINS: ANALYTIC 
AND NUMERICAL APPROACHES 

A. Analytic approach 

In this section, we consider line ar sp in chains with four 
spins. We follow Sect. IV of Ref. ^ (see also [U) and 
split the control problem for four spins into two subprob- 
lems for three spins (see Fig. 2]): The first subproblem is 
given on the first three spins by the time-optimal trans- 
fer from (1, 0, 0)'^ to (0, cos 7, sin 7)-^, where we are again 
using the coordinates (fi,r2,r3)^ of Sect. IIII Al Then, 
we apply certain (arbitrarily fast) hard pulses which can 
be easily determined by numerical methods. The second 
subproblem is given on the last three spins by the time- 
optimal transfer from (cos 7, sin 7, 0)"^ to (0,0,1)"^. In 
addition, we have to simultaneously search for the value 
of 7 which minimizes the pulse duration. This approach 
might not lead to time-optimal controls but simplifies the 
control problem significantly. 

The optimization of the considered subproblems can be 
reduced to time-optimal transfers from (cos a, sin a, 0)"^ 
to (0, cos/3, sin/3)^ for a,/? g [0,7r/2] generalizing the 
transfer of Sect. [IllAlfrom (1,0,0)'^ to (0,0,1)^. Us- 
ing methods of Ref. 12| we can find the optimal con- 




FIG. 6. The topology of the molecule "C^-^'^N-diethyl- 
(dimethylcarbonyl)fluoromethylphosphonate (see [13, [13) re- 
sults in coupling ratios fci = 2.38 and k2 — 0.94. Compare to 
Fig.H 



trols for the transfers using combined one-dimensional 
searches for the optimal initial values 9(0) and ^(0) 
of Eq. (g]). Both 6{0) and 9{0) are undetermined 
but related by 9{0) = sin[0(O)] cot a for the case of 
(ri(0), r2(0), 7-3(0))-^ = (cos a, sin a, 0)-^. The corre- 
sponding (semi-)analytic pulses are shown in Fig. [51 The 
values are motivated by the experimental system given 
in Fig. ini 



B. Numerical approach 

Motivated by the analytical approach, we numerically 
treat the control problem on four spins with two cases of 
coupling ratios (a) ki = 1 and fcj = 1 {J12 = J23 — 
J34 = 88.05 Hz), and (b) fci = 2.38 sa ^12/^23 and 
k2 ~ 0.94 PS J34/J23 (refer to Fig. [B] for the couphng 
values) . The coherence transfer is numerically optimized 
considering the following three levels of rf controls: First, 
we use only two different rf controls (one on each spin) 
operating on the second and third spin along the y- 
direction. Second, we use a total of four different rf 
controls (two on each spin) operating on the second and 
third spin along both the x- and j/-direction. Third, we 
use a total of eight different rf controls (two on each spin) 
operating on each of the four spins along both the x- and 
y-direction. The pulse shapes and the logarithmic TOP 
curves corresponding to two and eight rf controls (see Ta- 
ble IIII[) indicate that we do not gain a higher fidelity or a 
shorter duration from using more than the two controls 
f Table ITV)) . This is consistent with the analytic results, 
but the numerically-optimized pulses appear to be a little 
shorter than the analytical ones (cp. Fig. [S]). 



Conjecture 2 Consider a linear four- spin chain with lo- 
cal controls on each spin. One can time- optimally trans- 
fer coherence from Iix to SIiyl2yl5ylAz using only two 
controls along the y-direction which operate on the sec- 
ond and third spin, respectively. 



TABLE III. For linear four-spin chains with coupling ratios 
(a) ki — 1 and A;2 = 1 as well as (b) ki = 2.38 and ^2 — 0.94, 
the rf controls on the second and third spin along the j/-axis 
are plotted using solid-black and solid-gray lines, respectively. 
Other rf controls are plotted using a dashed-black line or in 
shades of gray. Compare to Table |I] 
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(b) ki = 2.38 and k2 = 0.94: 
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TABLE IV. For coherence transfers in linear four-spin chains 
(fci = fc2 = 1 as well as ki = 2.38 and ^2 = 0.94), we give 
the numerically-optimized times tp and the fidelities F in the 
case of two, four, and eight rf-controls. The duration tp is 
independent of the number of controls which suggests that 
only two rf-control on the second and third spin along the 
y-axis is sufficient for the time-optimal coherence transfer. 
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V. MORE GENERALLY-COUPLED SPIN 
SYSTEMS OF THREE AND FOUR SPINS 

In more generally-coupled spin systems, indirect cou- 
plings can strongly impede or enhance the coherence 
transfer. In this section we present detailed numerical 
optimizations and compare them to the case of linear 
spin chains. 

A. Three-spin system 

Along the lines of Section IIII Bl we numerically op- 
timize pulses for more generally-coupled three-spin sys- 
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tp (ms) 

conv. J,, = 88.05 Hz ... -conv. J.j = 70.00 Hz conv. /„ = 44.00 Hz conv. J,. = 0.00 Hz 

— .— oc. /„ = 88.05 Hz ■ oc. J, , = 70.00 Hz — •— oc. J,, = 44.00 Hz — •— oc. J,, = 0.00 Hz 

FIG. 7. (Color online) We compare numerically-optimized 
TOP curves (in shades of gray) for three-spin systems with 
one [Subfigure (a)] and six [Subfigure (b)] rf controls keeping 
Ji2 = J23 = 88.05 Hz constant while varying J13. At the same 
time we compare their performance with conventional pulse 
sequences (in shades of red) where both of the couplings J12 
and J13 evolve simultaneously. The limiting case of J13 = 
corresponds to the conventional pulse sequence of Sec. [Til Us- 
ing all six rf controls [Subfigure (b)] we can see higher fidelities 
F for smaller times and larger J13 compared to the case of 
only one rf control [Subfigure (a)]. A black arrow denotes 
where the numerically-optimized TOP curves for J13 — 44 Hz 
and J13 = Hz merge. 



tems keeping J12 = J23 ~ 88.05 Hz constant while vary- 
ing the additional coupling strength J13. By comparing 
the TOP curves for different values of J13, we conclude 
that for a larger coupling strength J13 the fidelity of the 
coherence transfer is smaller in the cases of one [Fig.[71[a)] 
and two (results are not shown) rf controls on the second 
spin. (As in Section llll Bl we obtain shorter pulse se- 
quences as compared to the conventional pulse sequence 
for Ji3 = 0.) However, using the rf controls on each 
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TABLE V. Numerical results for more generally-coupled 
three-spin systems with (a) couphng ratio k = 1 (i.e. J12 = 
J23 = 88.05 Hz) and additional couphng J13 — 2.935 Hz as 
well as (b) coupling ratio k = 1.59 (i.e. J12 = 73.1 Hz and 
J23 = 46 Hz) and additional coupling J13 = 10.0 Hz. Com- 
pare to Table m 



pulse shape 



logarithmic TOP curve 




(b) k = 1.59 

,60 
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of the three spins allows for a coherence transfer with 
higher fidelity while keeping the pulses short [Fig. [Tljb)]. 
Table |V] shows examples of shaped pulses and the corre- 
sponding logarithmic TOP curves for the coupling ratios 
k = 1 and k = 1.59. The coupling strengths are taken 
from the spin systems shown in Fig. [31 Detailed values 
are given in Table I VIII 



B. Four-spin system 



Following Section IIVBI we numerically optimize the 
shaped pulses for a more generally-coupled four-spin sys- 
tems. Analyzing the numerical results (see Table rvT| . 
we can say that this system needs all eight rf controls on 
each spin along both the x- and ^/-direction in order to 
achieve the coherence transfer with minimum duration 
and maximal fidelity. Table IVIII summarizes and com- 
pares the duration tp and fidelity F of shaped pulses for 
more generally-coupled spin systems. 



VI. EXPERIMENTAL RESULTS 

Analytic and numerically-optimized pulses are usually 
optimized for on-resonance cases. We follow the DANTE 
approach (60l . [61I |6^ in order to obtain pulses which are 
broadband, i.e. invariant with respect to the change of 



TABLE VI. Numerical results for a more generally-coupled 
four-spin systems with coupling ratios ki — 2.38 ~ J12/J23 
and k2 = 0.94 « J34/J23 (i.e. J12 = 46 Hz, J23 = 19.3 Hz, and 
J34 — 18.1 Hz) as well as additional couplings J13 =4.1 Hz 
and J24 = 2 Hz. Compare to Table iHll 



#u 



pulse shape 



logarithmic TOP curve 



U 200 
2 = -4001 



o -6 



16 36 
t (ms) 




57.8 20 45 

t, (ms) 
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70 



16 36 
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54 20 45 
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TABLE Vll. We compare the duration tp and fidelity F of 
numerically-optimized shaped pulses in the case of three- and 
four-spin systems allowing a varying number of rf controls 
u. Using only one rf control we show the effect of indirect 
couplings — which are usually present in experiments — on the 
fidelity of optimized pulses. Hence more rf controls are nec- 
essary for higher fidelities. The J-values are taken from the 
actual spin systems shown in Figs. O and [S] 
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the chemical shift in a given offset range (Fig. |8]). First, 
a shaped pulse is converted into a sequence of short hard 
pulses and delays. We used hard pulses with constant 
flip angles (see below), and the delays between the hard 
pulses correspond to the time required by the shaped 
pulse to accumulate this flip angle [60l. [61I [63 | . Then, a 
rcfocusing element (i.e. 7r-pulse) j29l. l6l| is inserted be- 
tween two hard pulses. The offset bandwidth covered by 
a refocused DANTE sequence is directly proportional to 
the rf amplitude of the hard pulses and the 7r-pulses used 
in the sequence. 

All the experiments are implemented on a Bruker 
AVANCE III 600 MHz spectrometer at 298 Kelvin: We 
use a triple resonance TXI probe head with Z-gradient 
in the case of the three-spin system with fc = 1. For the 



r (a) 



A, 


A, 


4 




A, 


A„ 


A, 1 A, 


A.| A. 


A. 1 A, 


A, 1 A 


A, 1 A 


A. 1 A, 



FIG. 8. In the DANTE approach an on-resonance shaped 
pulse (see (a)) is converted to a series of short hard pulses 
and delays Ai (see (b)). Then, the pulse can be converted 
to a broadband pulse by inserting a refocusing element (i.e. 
TT-pulse) represented by solid bars between two hard pulses 
(see (c)). 
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Experiment (a) Simulation (a') 




Av (kHz) Av (kHz) 



FIG. 10. (Color online) We compare for a three-spin system 
the offset (Av) profile for ±2 kHz of the anti-phase signal 
(see text) resulting from a conventional pulse sequence {tp = 
17.7 ms) in the case of experiment (see (a)) and simulation 
(see (a')) with broadband versions of the analytic pulses (tp = 
15.5 ms, see (b) and (b')) for the case of coupling ratio k = 
1.59. 



Experiment (^a) Simulation (a') 




-1 1-1 1 

Av (kHz) Av (kHz) 



FIG. 9. (Color online) We compare for a three-spin system 
the offset (Av) profile for ±1 kHz of the anti-phase signal 
(see text) resulting from a conventional pulse sequence {tp = 
11.4 ms) in the case of experiment (see (a)) and simulation 
(see (a')) with broadband versions of the analytic pulses {tp — 
9.8 ms, see (b) and (b')) for the case of coupling ratio k = 1 
of three spins. 



three spin system with k = 1.59 and the four spin sys- 
tem with ki = 2.38 and k2 — 0.94 we use a custom-made 
6-channel probe head with Z-gradient addressing all nu- 
clei i^F, iH, 3ip, 12c (or 13C), and (or ^^N) (see 
[13, nil)- In the experiments for three spins we use the 
molecules shown in Fig. [3] The experiment for the cou- 
pling ratio k = 1 uses the first molecule (see Fig. ^a.)) 
which is dissolved in deuterated water D2O. For k = 1.59 
we use the second molecule (see Fig. Elb)) dissolved in 
deuterated methanol CD3OD. The simulated and ex- 
perimental offset profiles are shown in Figs. [S] and [TUl 
We emphasize that the duration of the broadband ver- 
sions of the analytic or the numerically-optimized pulses 
is shorter than for the conventional pulse sequence while 
keeping its robustness. 

We first discuss two three-spin systems: In the case 
of fc = 1, we start from the initial polarization /12 of 
(which models the first spin) and apply a ^-pulse 
along -|-?/-direction in order to obtain the coherence Ii^- 
By applying a broadband version of our shaped pulse 
of ^^N (which models the second spin) we get the three- 



spin coherence 'iliyl2yl3z- The broadband version of this 
shaped pulse is divided into four hard pulses with an 
amplitude of 4145.936 Hz, a flip angle of 45.00 degrees, 
and zero phase; it also contains refocusing 7r-pulses where 
the phases are chosen according the MLEV-4 cycle [65j . 
Next, we apply ^-pulse on ^^N along the x-direction and 
we obtain the coherence ^Iiyl2zl3z- In the end, we can 
detect an anti-phase signal of -'^H (first spin) with respect 
to the spins of ^^N and (which models the third spin). 
Similarly, in the case of fc = 1.59, we start from the initial 
coherence Iiz on the spin of (which models the first 
spin) and apply a -l-pulse along -t-y-direction in order to 
obtain the coherence Iix- Then, we apply the broadband 
version of our shaped pulse on the spin of ^^F in order to 
produce the three-spin coherence 'iliyl2yhz- The broad- 
band version of this shaped pulse is divided into four 
hard pulses with an amplitude of 10000 Hz, a fiip angle 
of 45.03 degrees, and zero phase; it also contains refo- 
cusing TT-pulses where the phases are chosen according 
the MLEV-4 cycle [g^. In the next step, we apply a ^- 
pulse on the spin of ^H and we end up with the coherence 
4/12/2^/32. Finally, we detect an anti-phase signal on the 
spin of ^^F with respect to the spins of ^H and '^^P. 

In the four-spin system, we show on-resonance simula- 
tions and experiments for numerically-optimized shaped 
pulses comparing the conventional approach with ana- 
lytical and numerically-optimized pulses (see Fig. [TT]) . 
The corresponding experiments are implemented on the 
molecule of Fig. [51 which is dissolved in deuterated ace- 
tonitrile. FiguresfTTTb) and (c) show a reduction in signal 
intensity for the simulation if we compare the effect of the 
pulse on the abstract linear spin chain (shown in gray) 
with the effect on the more realistic and more generally- 
coupled spin system (shown in black) as the correspond- 
ing pulses were only optimized for the abstract linear 
spin chain. We remark that the pulse of Fig. rTTT d) is op- 
timized for a more generally-coupled spin system while 
using rf-controls on all spins. Thus, we conclude — using 
also the data of Table Elh-that the pulse of Fig. WVid) 
shows a higher fidelity when compared to the pulses of 
Figs. fTlTb) and (c). Furthermore, the pulse correspond- 
ing to Fig. [TTT d) is shorter (by 14%) than the convcn- 
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FIG. 11. We show the anti-phase signal of the spin of ^'"'F 
with respect to the spins of ^H, ^^C, and ^^N of a four-spin 
system corresponding to the simulation (left) and experiment 
(right). We use the conventional pulse sequence (a) {tp — 
64.4 ms), an analytical pulse sequence (b) {tp = 53.9 ms), a 
pulse which was numerically-optimized for the abstract linear 
spin chain with rf controls on the second and third spin (c) 
{tp — 53.2 ms), and a pulse which was numerically-optimized 
for the more generally-coupled spin system with rf controls on 
all spins (d) {tp = 54.0 ms). The simulation for the abstract 
linear spin chain is given in gray color. All the other plots for 
the more realistic case of a more generally-coupled spin system 
are given in black color. The plots are scaled vertically by a 
factor of two. 



tional pulse sequence corresponding to Fig. [TlTa) while 
maintaining its robustness to additional couplings (see 
also Tabic EmJ. 



VII. LINEAR SPIN CHAINS WITH MORE 
THAN FOUR SPINS 



In this section, wc generalize the numerical optimiza- 
tion of shaped pulses to linear spin chains of five and more 
spins. Figure [T^] shows two examples of the optimized 
pulse shapes with coupling ratios k{ = 1 and coupling 
strengths J^^e+i = 88.05 Hz. These examples suggest 
that time-optimal controls can be obtained on multiple 
spins even while irradiating only on the spins two to £—1 
along the y-direction (cp. Sec. IV of Ref. [T^l)- We obtain 
shorter pulses for the numerically-optimized pulses com- 
pared to the conventional pulse sequences as summarized 
in Table Eini 




6 12 
t (ms) 



21.06 



FIG. 12. Using all ten (or twelve) rf controls we determined 
numerically-optimized pulse shapes for linear spin chains of 
length five (or six) in the case of fc = 1. We remark that most 
control strengths are very small. 



TABLE VIII. We compare the minimum time tp required for a 
coherence transfer by numerically-optimized (oc) and conven- 
tional (conv) pulse sequences for different number ?i of spins 
and coupling ratios k. 
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^ For n = 3 we have k - 
and k2 = 0.94. 



1.59. And for n = 4 we have fci = 2.38 



VIII. CONCLUSION 

In the case of linear three-spin chains we reproduced 
numerically the previous analytical results [III, ob- 
taining the same family of restricted controls by applying 
pulses only on the second spin along the y-axis. The same 
holds for linear four-spin chains where we also obtain the 
analytical family of restricted controls by applying pulses 
only on the second and third spin along the y-axis; but 
the numerically-optimized pulses appear to be a little 
shorter than the analytical ones. Both for three and four 
spins no gain in pulse duration is found if arbitrary pulse 
structures are allowed. These observations are summa- 
rized in Conjectures [Hand [2] Even for longer spin chains 
(consisting of up to six coupled spin-1/2) there is some 
numerical evidence suggesting that the same restricted 
controls motivated by Refs. [ll|, [T^] lead to time-optimal 
pulses (in the case of unrestricted controls) for linear spin 
chains of arbitrary length. 

Further numerical results are presented for more gen- 
eral and more realistic coupling topologies, for which so 
far no analytical results are known. Compared to linear 
spin chains we obtain different pulse structures depend- 
ing on the number of available controls. We hope that the 
presented results and conjectures will motivate further 
analytical work in order to develop a better understand- 
ing of time-optimal control sequences for the generation 
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of multi-spin coherence. 

Note that the minimum times for the transfers lis ~^ 
2^~^Iiei ■ ■ ■ Ine„ 8-11(1 2"~^/iej^ • • ■ /„g^ lis are identical 
((5, efc € {x, y,z} for k = 1, . . . , n), which is directly rel- 
evant for "out and back" experiments and the reconver- 
sion of multiple-quantum coherence to detectable single 
quantum operators. In the experimental part, we demon- 
strated that the optimized pulse sequences work in real- 
istic settings under relaxation and experimental imper- 
fections (e.g. inhomogeneity of the control field, miscal- 
ibrations, and phase transients). In addition, the pulses 
can be made broadband (i.e. robust with respect to fre- 
quency offsets) using the DANTE approach. 

Here we assumed for simplicity that each spin- 1/2 
can be selectively addressed, which is directly relevant 
to heteronuclear spin systems but the optimal transfer 
scheme can also be adapted to homonuclear spin sys- 
tems. The presented sequences can be directly applied 
to small molecules and peptides, which is in particular 
true for the broadband versions. The minimum pulse 
sequence durations for complete transfer are reduced by 
up to 24% compared to conventional approaches (see Ta- 
ble IVIIip . Conversely, for a fixed transfer time signif- 



icantly improved transfer amplitudes are possible, e.g., 
for a linear three-spin chain we gain approximately 23% 
in transfer efficiency when we allow only for half of 
the transfer time necessary for a complete transfer (cp. 
Fig. 7). For large proteins, further gains in efficiency are 
expected if relaxation-optimized pulse sequences can be 
developed for the specific relaxation super operator given 
in the system. Although such sequences are beyond the 
scope of the present paper, the results on time-optimal se- 
quences presented here provide an important benchmark 
for relaxation-optimized sequences. 
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